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Abstract 

We present the results of exact diagonalization calculations of the isolated and isothermal on-site static susceptibilities in the 
anisotropic extended Heisenberg model on a linear chain with periodic boundary conditions. Based on the ergodicity considerations 
we conclude that the isothermal susceptibility will diverge as T — > both in finite clusters and in the bulk system in two non-ergodic 
regions of the phase diagram of the system. 
© 2008 Elsevier B.V. All rights reserved. 
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It is often thought that the ergodicity - the property to 
explore the whole Hilbert space during time evolution - 
is a natural attribute of a physical system and, moreover, 
that its possible violation is hard to observe. Ergodicity as- 
sumptions are usually made when passing through various 
self-consistent schemes in determining susceptibilities, cor- 
relation functions etc. Nevertheless, non-ergodic behavior 
may manifest under certain circumstances. The purpose of 
this article is to show how non ergodicity can significantly 
change the properties of the system and may become in 
principle observable. 

When dealing with the static response of a quantum sys- 
tem two definitions of the susceptibility are usually con- 
sidered: the static "isolated", or "Kubo susceptibility", %o 
and the isothermal one \ T [!]• These two definitions cor- 
respond to different environmental conditions under which 
the response is measured. Namely, xo is derived upon the 
assumption that the system is isolated, while in order to 
measure x T one nas to maintain the system in thermal 
equilibrium with an external bath at a given temperature. 
As it was pointed out by Kubo and Suzuki [2] , a necessary 
and sufficient condition for these two response functions to 
coincide is the ergodicity of the system. By considering as 
a specific example the response function of an operator A 
to a perturbation coupled to A itself, the isolated Kubo 
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susceptibility will be: 

Xo = - i 
while the isothermal one: 



dt([A(t),A(0)}) 



(i 



r = I d\(A(-i\)A)-p(Ay 

la 



Their difference, as shown in [1,2], is given by: 

x T - xo = /?(r - r«*), 

where we have introduced the quantity T: 



T = lim (A(t)A(0)) 

t — >-oo 



and its ergodic value 
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When r ^ T erg the response is said to be non-ergodic and 
the two susceptibilities differ. In this case, non-ergodicity 
becomes observable as it provides a diverging contribution 
to x T in the limit T -> 0. 

In this paper we show the results for the isolated suscep- 
tibility in the extended anisotropic Heisenberg model on a 
linear chain of 14 sites with periodic boundary conditions. 
The susceptibility has been calculated using the spectrum 
of the system obtained by means of the exact numerical 
diagonalization. The Hamiltonian of the system reads as: 
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Fig. 1. On-site magnetic susceptibility in the various regions from the phase diagram (see [4]) as a function of temperature for the 
Hamiltonian (6) on a linear chain of 14 sites. Panel a) collects the results for the points from Ergodic I (solid and dotted lines) and Ergodic 
II phases (dashed and dashed-dotted lines) while the panel b) those from Non-Ergodic I (solid and dotted lines) and the Non-Ergodic II 
phases (dashed and dashed-dotted lines). 
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We fix J z > so that the z-component of the nearest- 
neighbors exchange is ferromagnetic and use J z as the en- 
ergy unit. We use S" as the operator A in (1): 



Xo(M) = -i/ dt([Sf(t),S?(0)}). 



(7) 



Because of the translational invariance Xo(h i) is indepen- 
dent of the site i. 

Previously [3], [4] we have already studied the question 
of ergodicity of the response of S 1 ? in the Hamiltonian (6) . 
We have constructed a phase diagram of (6) in the plane 
J' — J± for both zero and finite temperatures using as large 
as 26 sites clusters. We have concluded that at T = and 
in the bulk limit there are two non-ergodic regions and two 
ergodic ones (see Fig.l in [4]). A transition zone exists at 
finite sizes, which probably would become a transition line 
in the bulk system. At T > the finite-size scaling indicated 
that the system is always ergodic. 

On Fig. 1 the Kubo susceptibility is plotted as a function 
of temperature for the ergodic (panel a) and non-ergodic 
(panel b) regions from the phase diagram obtained in [4] . 
While at high temperatures Xo(T) ~ 1/T, as it should 
be after the Curie law, the most interesting behavior is 
concentrated at low temperatures. One can easily show that 
when T — > the Kubo susceptibility can be rewritten as: 



Xo(T) = 



9 N 

-Y Y 

1=1 E n > E 



\(0,l\Sf\n}\< 
E n — E 



+ 0(T), (8) 



where N is the number of degenerate ground states (|0, 1), 
1 = 1,... N), \n) is the n— th excited eigenstate of (6) and 
f2(T) goes exponentially to zero when T — > 0. (8) is nothing 
else but the result of the perturbation theory up to the sec- 
ond order. One can see from (8) that Xo(0) is non-singular 
and in general non-zero since all the terms under the sum 



have the same sign. This is not the case in the phase Non- 
ergodic I (see the dashed and dashed-dotted lines on the 
Fig. lb)) since in this phase the ground state is doubly de- 
generate with all spins either up or down. Such states are 
eigenstates of Sf and therefore the matrix elements in (8) 
vanish. All the other phases have their ground states con- 
nected to the rest of the Hilbert space by 5? so that for 
them xo(0) # 0. 

From our data for the isolated susceptibility and by us- 
ing (3), we can find the temperature dependence of \ T as 
well. It is clear from (3) that in both ergodic phases (Er- 
godic I and Ergodic II) x T {T) = Xo{T). On the contrary, 
the isothermal susceptibility will diverge at low T as T/T, 
where T = 1/12 + 1/61/ for the phase Ergodic I on a cluster 
of L sites, and T = 1/4 for the phase Ergodic II, indepen- 
dently on L. Based on the finite-size scaling for T made in [4] 
these results will hold also in the thermodynamic limit. 

In conclusion, we have studied the on-site static magnetic 
response of (6), considering both definitions of the suscep- 
tibility. While for the isolated one the non ergodicity enters 
implicitly through the degeneracy of the eigenvalues, the 
isothermal susceptibility contains an explicit term diverg- 
ing at low-T in the case of non-ergodic phases. In real mate- 
rials, however, this divergence might not be observable due 
to perturbations such as dipole-dipole interactions which 
destroy the conservation of z-component of the total spin 
and ensure the ergodicity. 
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